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O . We present Suzaku observations of five hard X-ray selected nearby Seyfert 2 

galaxies. All the sources were clearly detected with the pin Hard X-ray Detector 
up to several tens of keV, allowing for a fairly good characterization of the broad- 
band X-ray continuum. We find that a unique model, even including multiple 
components, fails to represent the spectra of all the sources. Heavy obscuration 
manifests itself in different flavours. For two sources there is evidence for a 
reflection dominated continuum; among the other three, one is "mildly" Compton 
thick [Nh ~ 10^^ cm~^), while the remaining two are heavily obscured {Nh — 
]^g23.5 cj-Q-2~)^ ]^-^J^^ Comptou thin. Strong, narrow, iron Ka lines {EW ~ 1-2 keV) 
due to neutral or mildly ionized gas, are detected in Compton thick AGN. In 
all of them the Ka line is accompanied by the K/3. The intensity and shape of 
the soft X-ray spectrum are different from object to object. Soft X-rays may 
originate from a nuclear component scattered off, or leaking through, the X-ray 
absorber, plus thermal X-rays from the host galaxy. Emission from circumnuclear 
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gas photoionized by the active nucleus, parameterized with a power law plus 
individual narrow Gaussian lines, also provides an acceptable description of the 
soft X-ray spectra. The limited Suzaku XIS CCD energy resolution does not allow 
us to draw firm conclusions on the origin of the soft X-ray emission. We briefly 
discuss our findings in the light of AGN Unified model and the geometry of the 
obscuring gas. 

Subject headings: galaxies: active — X-rays: observations — Seyfert galaxies: 
individual (NGC 5728, NGC 4992, ES0263-G13, ESQ 137-G34, ES0323-G32) 



1. Introduction 

A fraction as high as 50% of Seyfert 2 galaxies in the nearby Universe are obscured 
in the X-ray band by column densities of the order of, or larger than the inverse of the 
Thomson cross-section [Nu > a^^ :^ 1.5 x 10^'* cm~^), hence dubbed Compton thick (CT). 
If the optical depth (r = Nhctt) for Compton scattering does not exceed values of the 
order of "a fetif\ X-ray photons with energies higher than 10-15 keV are able to penetrate 
the obscuring material and reach the observer. For higher values of r, the entire X-ray 
spectrum is depressed by Compton down scattering and the X-ray photons are effectively 
trapped by the obscuring material irrespective of their energy. The former class of sources 
(mildly CT) can be efficiently detected by X-ray instruments sensitive above 10 keV, while 
for the latter (heavily CT) their nature may be inferred through indirect arguments, such 
as the presence of a strong iron Ka line over a fiat reflected continuum. The search for 
and the characterization of the physical properties of CT AGN is relevant to understand 
the evolution of accreting Supermassive Black Holes (SMBHs). In particular, mildly CT 
AGN are the most promising candidates to explain the so far largely unresolved spectrum of 
the X-ray background around its 30 keV peak (Worsley et al. 2005; Treister & Urry 2005; 
Gilli et al. 2007). According to the Gilli et al. (2007) XRB synthesis model, their integral 
contribution to the hard X-ray background is of the order of 25-30%. This fraction has 
been estimated under simplified hypotheses. In particular, the same luminosity function and 
cosmological evolution of unobscured and Compton thin AGN is assumed; moreover, the 
number density of "heavily" CT AGN is the same of "mildly" CT. While these assumptions 
are not inconsistent with the present observational framework (Risaliti et al. 1999; Guainazzi 
et al. 2005), it should be noted that absorption column densities in excess of 10^^ cm~^ were 
measured or inferred for about a few tens of nearby AGN (Comastri 2004; Delia Ceca et 
al. 2008), and only a handful of them are known beyond the local Universe (Norman et al. 
2002; Iwasawa et al. 2005). CT AGN may represent a significant fraction of the accretion 
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power in the Universe and indeed a correction for the CT contribution to the estimates of 
the SMBH local mass density has to be included in the calculations (Marconi et al. 2004). 

An unbiased census of extremely obscured AGN would require to survey the hard X-ray 
sky above 10 keV with good sensitivity. Such an argument is one of the key scientific drivers 
of the NUSTAR (Harrison et al. 2010) and astro-h (Takahashi et al. 2010) missions, which 
will be launched in the next few years, and of the nhxm mission study (Pareschi et al. 2010). 

For the time being, one has to rely on the observations obtained by the high-energy de- 
tectors on board BEPPOSAX, INTEGRAL, SWIFT and, more recently, Suzaku. All-sky surveys 
were performed using both the IBIS coded-mask telescope onboard INTEGRAL and the BAT 
detector on board SWIFT. Though limited to bright and thus low-redshift sources, they have 
proven to be quite successful. More than one hundred AGN are reported in both integral 
(Beckmann et al. 2009) and swift catalogues (Tueller et al. 2010; Cusumano et al 2010). 

Hard X-ray selection is less biased against absorption and thus provides a useful bench- 
mark to study the column density distribution of obscured AGN and the fraction of CT 
sources in the local Universe. The large majority of INTEGRAL and SwiFT sources were 
observed by XMM and Chandra and their broad band X-ray spectra discussed by Winter 
et al. (2008). While the already known, nearby CT AGN were recovered by Swift and 
Integral surveys, the fraction of "newly" discovered CT AGN is surprisingly low and ap- 
parently inconsistent, by about a factor of 2, with that predicted by Gilli et al. (2007) in the 
local Universe. On the basis of these findings, it has been proposed (Treister et al. 2009) 
that the contribution of CT AGN to the hard X-ray background may be significantly lower 
(by a factor of 2 to 3) than previously thought, with important implications for the evolution 
of the accretion power. 

Since absorption column densities are often measured on relatively poor quality X- 
ray spectra taken from non-simultaneous observations, above and below about 10 keV, and 
combining different instruments, the CT fraction and, more in general, the absorption dis- 
tribution in the local Universe may still be subject to several uncertainties. Good quality 
simultaneous spectra extending over the 0.5-100 keV energy range are needed for a robust 
measurement of absorption column densities, especially in the Compton thick regime. More- 
over, nearby obscured AGN always show excess emission above the extrapolation of the 
obscured nuclear spectrum (e.g. Turner et al. 1997). Soft X-ray spectroscopy is a powerful 
tool to study the origin of this component which is plausibly related to warm gas photoion- 
ized by the nuclear continuum (Guainazzi & Bianchi 2007). The X-ray detectors onboard 
the Japanese satellite Suzaku are well suited to this purpose, at least as far as "local" AGN 
are concerned. We have conceived a program with Suzaku to observe nearby, relatively X-ray 
bright (> 10~^^ erg cm~^ s~^) AGN selected by requiring i) a significant detection above 10 
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keV by either Swift or Integral and ii) evidence for significant X-ray absorption at lower 
energies from archival observations. The immediate science objective is the characterization 
of the broad band X-ray spectra of hard X-ray selected, heavily obscured AGN and a first 
step towards a better census of Compton thick absorption in the nearby Universe. 



2. Sample selection 

Five sources (NGC 4992, NGC 5728, ES0137-G34, ES0263-G13, ES0323-G37) spec- 
troscopically classified as Seyfert 1.9-2.0 galaxies, were originally selected from the inte- 
GRAl/ibis (Beckmann et al. 2006) and swift/bat (Markwardt et al. 2005) catalogues. 
The column densities, as inferred from archival ( Chandra and XMM-A/^ewton) observations, 
are of the order of 10^^~^*^ cm~^, though affected by large errors. By construction, the sample 
is biased towards obscured AGN. 

The extremely bright flux of NGC 4992 in the INTEGRAL observation, coupled with the 
lack of a detection in the RASS (Voges et al. 1999), is consistent with an absorbed spectrum, 
as confirmed by a snapshot Chandra observation (Sazonov et al. 2005). 

Chandra revealed an heavily obscured nucleus in NGC 5728 (Zheng et al. 2006). Al- 
though the column density is not well constrained due to the limited response at high energies, 
the presence of a strong iron line (EW ~ 1 keV) suggests that the source might be obscured 
by Compton thick gas. 

ES0137-G34 is detected by Integral and subsequently observed with IQAVi-Newton 
(Malizia et al. 2009). The combined XMM and Integral spectrum can be well fitted 
by both a transmission and a reflection model as well as by a more complex absorption 
distribution. In all cases, there is evidence for Compton thick absorption. 

ES0323-G32 and ES0263-G13 were selected on the basis of a significant (> 5(t) inte- 
GRAl/ibis detection, unambiguous optical classification as Seyfert 2 galaxies, and the lack 
of archival observations in the 2-10 keV energy range. Both are detected by rosat with 
a soft X-ray spectrum. The extrapolation of the best-fit integral spectrum in the soft 
ROSAT band suggests obscuration in excess of 10^'^ cm~^. 



3. Observations and data Analysis 

The fifth Japanese X-ray satellite Suzaku (Mitsuda et al. 2007) carries four X-ray 
telescopes with a CCD camera, the X-ray imaging spectrometer (XIS; Koyama et al. 2007) 



-5- 



in the focal plane of each of them, and a non-imaging Hard X-ray Detector (HXD; Takahashi 
et al. 2007), consisting of Silicon pin photo-diodes and GSO scintillation counters. They 
simultaneously cover the 0.2-12 keV and 10-700 keV energy ranges, respectively, although 
with different sensitivities. 

The XIS cameras were operating with the Space-Row Charge Injection in the standard 
mode. The data were taken by the three XISs (two FI CCD cameras - XISO and XIS3 - 
and one BI CCD camera XISl) and the HXD pin. For the two sources - NGC 5728 and 
NGC 4992 - observed at the beginning of the mission, data from the FI CCD XIS2 are also 
included. All the sources were observed in the "HXD" nominal position on the detector 
plane which is optimised for the maximum throughput of the HXD rather than the XIS. 



The data from version 2 processing^ were used for the analysis. The cleaned event files 
produced by the data processing with the standard selection criteria were used. The non 
X-ray background of the pin detector was estimated by the "tuned" background model. 
The effective exposure times and the background subtracted count rates in the xis and pin 
detectors for each object are reported in Table 1. 

Spectral fits are performed with XSPEC version 11.3. Front Illuminated CCD were 
co-added, while the back illuminated chip was discarded after verifying that its spectrum 
does not significantly improve the quality of spectral fits. 

The data from XIS and pin detectors were fitted together. The intercalibration constant 
(pin/XIS) is l.lqj. The broad band Suzakuspectra. in a vF^ representation of the five sources 
are shown in Figure 1. In Figure 2 we report the XIS spectra. 



4. Spectral analysis 

For sake of clarity and to ease the comparison with some of the previous Suzaku results 
on heavily obscured AGN (Ueda et al. 2007; Eguchi et al. 2009), we consider two baseline 
spectral models for the broad band Suzaku spectra, namely a transmission dominated (TD) 
model: waca^Capec + pexrav + zpcafbs (cutoff pl+2xzgauss) ) in XSPEC notation and 
a reflection dominated (RD) model: waca^Capec + pow + pexrav + 2xzgauss). The spec- 
tral components of the TD model are: a power law with an exponential cut-off fixed at 200 
keV (Gilli et al. 2007) for the primary spectrum plus intrinsic absorption, a reflection com- 



^http:// www.astro.isas.ac.jp/suzaku/process/caveats 

^ftp:// legacy.gsfc.nasa.gov/suzaku/doc/xrt/ suzakuniemo-2008-06.pdf 
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ponent from neutral gas (pexrav; Magdziarz & Zdziarski 1995) and iron K emission (both 
Ka and K/3, whenever required by the spectral fits). The photon index of the reflection 
component is linked to that of the primary continuum and no absorption on the reflected 
spectrum is considered. The RD model includes the same spectral components except for 
the absorbed primary power law. The photon index of the RD model is left free to vary. 
The intensity of the reflection component R in TD model is defined as the ratio between 
the normalization of the reflected and primary (i.e. power law) spectrum. The inclination 
angle of the reflection component is fixed to cosi=0.9, which is equivalent to assume a rather 
face-on geometry for the reflecting material, likely to be associated with the inner walls of 
the obscuring torus. 

Soft X-ray emission in excess of the extrapolation of the absorbed nuclear continuum 
at low energy is observed in all the sources of our sample. This emission may be due to hot 
gas associated to the host galaxy and is modeled with an updated version of the original 
Raymond & Smith (1977) code for thermal plasmas (apec in Xspec), and to AGN primary 
continuum scattered on the line of sight or leaking through the nuclear absorber, represented 
by a partial covering (zpcfabs in XSPEC), in the TD model, and with a power law in the 
RD model. In both models the power law slope is assumed to be the same of the intrinsic 
hard X-ray continuum, with the exception of one object (see § 4.4). The scattering fraction 
/ can be defined only for the TD model, because the level of the intrinsic continuum in RD 
fits is unknown. In order to compare the intensity of the soft X-ray power law in TD and 
RD models, we report the monochoromatic flux densities at 1 keV in Table 2. 

We have also considered an alternative possibility for the origin of soft X-ray emission 
in obscured AGN in terms of emission from a photoionized plasma modeled in XSPEC as 
wacrajCpow + Nxzgauss + wa(cutoffpl + 2xzgauss)) for TD sources and as wacJa/Cpow 
+ Nxzgauss + pexrav + 2xzgauss)) for RD AGN (Table 3). A^ is the number of fitted 
lines in the soft X-ray band which varies form source to source (Table 4). The additional 
zgauss component accounts for Ka and K/3 lines. Given that the XIS CCD energy resolution 
is not appropriate to fit a photoionization code model (Iwasawa et al. 2003), the expected 
soft X-ray spectrum of a photoionized plasma is modeled with a power law with photon index 
free to vary (Table 3), plus intrinsically narrow (setting a=0 in the zgauss XSPEC model) 
Gaussian lines at the energy expected (Table 4) for the strongest features observed in high 
resolution RGS spectra of obscured AGN (see Guainazzi & Bianchi 2007). The presence of 
the line is evaluated on the basis of a visual inspection of the residuals obtained by fitting a 
single power law. The best fit energies, fluxes and associated la errors of the various lines 
for each source are reported in Table 4. 

The best fit spectral parameters of the continuum along with observed X-ray fluxes and 
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2-10 keV intrinsic luminosities are summarized in Table 2 for the fits with a partial covering 
and a thermal spectrum for the soft X-rays, and in Tables 3 and 4 when the soft X-rays are 
modeled with a power law plus Gaussian lines. The shape of the hard X-ray continuum is 
rather independent on the model adopted for the soft X-ray spectrum. Therefore fluxes and 
luminosities are reported only in Table 2. Galactic absorption using Morrison & McCammon 
(1983) cross sections (model vsicai in XSPEC) is fixed at the values measured by Dickey & 
Lockman (1990). 

The detailed results of spectral fits for each source are summarized in the following. 



4.1. NGC 4992 

The TD model provides an acceptable fit to the observed spectrum of NGC 4992, while 
the RD model is ruled out at an extremely high confidence level (Ax^ > 600). The primary 
continuum is obscured by cold gas, with a column density of N^ ~ 5.5 x 10^^ cm~^ almost 
fully covering the central source. The iron line intensity {EW ~ 340 eV) is consistent with 
being originated by transmission through the obscuring gas. The constraints on the fraction 
of the scattered {fscatt < 0.3 %) and refiection intensity (R ~ 0.2-0.5) are rather stringent. 
The soft X-ray emission is extremely weak, and both a thermal spectrum {kT ~ 0.3 keV) and 
a steep power law plus three emission lines at the energies corresponding to the transitions 
reported in Table 4 provide an acceptable fit. We note that in the latter case no continuum is 
formally required, and the entire soft X-ray fiux may be explained as emission from multiple 
lines. 



4.2. NGC 5728 

The TD model provides the best fit to the broad band spectrum of NGC 5728. A RD 
model can be safely discarded on a statistical basis (A^^ > 130). The source is mildly 
Compton thick {Nh c^ 1.4 x 10^^ cm^^) and thus hard X-rays (> 10 keV) are piercing 
through the obscuring gas. The reflection component is rather weak {R c^ 0.1-0.3). A 
strong (EW ~ 1 keV) iron Ka line is detected along with a line at ~ 7 keV. The best fit 
energy of the ~ 7 keV feature, its intensity (EW ~ 70-90 eV) and the observed intensity 
ratio with respect to Ka, favour an interpretation in terms of K/3 emission. The soft X-ray 
spectrum below a few keV is fitted by a combination of a scattered power law, with fscatt — 
0.4-1.2 % , and thermal plasma emission {kT ~ 0.3 keV). The residuals below a few keV 
suggest the presence of line emission. A fit with a steep power law plus several emission lines 



(Table 4) also provides a statistically acceptable solution (Table 3) and smooth residuals in 
the soft X-ray band. 



4.3. ESO 263-G13 

The Suzaku spectrum of ES0263-G13 is best fitted by a TD model. A RD model fails 
to fit the data over the entire energy range (A^^ > 700). The absorption column density is 
the lowest in the present sample (~ 3 x 10^^ cm~^); the scattering fraction and the refiection 
component intensity are poorly constrained {fscatt < 0.8 % and R ^ 0.1-1.2, respectively). 
The iron Ka line EW (~ 80 eV) is consistent with that expected by transmission through 
the observed column density. The soft X-ray spectrum is characterized by low counting 
statistics. Both a power law plus thermal X-rays {kT ~ 0.5 keV) and a power law plus 
two emission lines provide an acceptable fit. The latter fit is slightly better in terms of x^, 
though the improvement is not statistically significant. 



4.4. ESO 137-G34 

Both a RD model {x^/d.o.f ~ 169/139 see Table 3) and a TD model {x^/d.o.f ~ 
160/138) provide an acceptable fit to the broad band spectrum of ESO 137-G34. In the 
TD fit, the primary hard X-ray continuum (F ~ 1-8), is obscured by a column density of 
~ 1.2 X 10^^ cm^^. The intensity of the scattering fraction {fscatt — 5%) is higher than 
that measured in the sources best fitted with a TD model, while the refiection component is 
relatively weak {R ~ 0.3). Even though none of the two models cannot be preferred on the 
basis of purely statistical arguments, the individual parameters of the TD model are highly 
degenerate, and cannot be reliably constrained. As a consequence, a RD model is preferred. 
It is worth noting that the source would be classified as Compton Thick in both cases. The 
strong Ka line (EW ~ 1.5 keV) is accompanied by a line at ~ 7 keV. The best-fit line 
intensity ratio (~ 0.2) is slightly higher, but not inconsistent within the errors, than that 
expected by a K/3 line origin for the ~ 7 keV feature. The soft X-ray spectrum can be fitted 
with a thermal plasma {kT ~ 0.8 keV) and a power law slope which is significantly steeper 
than that of the hard X-ray continuum. Some residual line like emission, around 1 keV and 
1.8 keV, is still present, suggesting a more complex spectrum. A fit with a power law plus 
several Gaussian lines (Table 4) provides an acceptable solution in terms of x^ statistic and 
smooth residuals (see Fig. 4). 
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4.5. ESO 323-G32 

Both a RD model {x^/d.o.f ~ 79/89 see Table 3) and a TD model {x^/d.o.f ~ 83/88) 
provide an acceptable fit to the broad band spectrum. In the TD fit, the primary hard X-ray 
continuum F ~ 2.0, is obscured by a column density of the order of ~ 1 — 2 x 10^^ cni"^. 
The relative intensity of the transmitted component is about a factor two lower than that 
of the reflection component. Similarly to ESO 137-G34, the best fit spectral parameters, in 
particular the column density Nh and the scattering fraction, cannot be constrained. A RD 
model is thus preferred as the best fit to ESO 323-G32, the faintest source of the sample. 
The iron Ka line is very strong (EW ~ 2 keV), and the best fit energy and intensity of 
the ~ 7 keV feature are consistent with a K/3 origin. The weak soft X-ray emission can be 
adequately modeled with either a power law plus a thermal spectrum or with a power law 
plus two Gaussian lines. Both fits are statistically acceptable and comparable. 



4.6. Summary of fitting results 

4.6.1. Hard X-ray spectra 

All the targets are obscured by column densities larger than ~ 3 x 10^^ cm^^. Two of 
them are highly absorbed but Compton thin (A^^^;^ ~ 3 — 5 x 10^^ cm~^), while the other 
three are Compton Thick. In NGC 5728 the primary hard X-ray continuum is transmitted 
through a column density Nh ^ lO^^ cm'^. In ESO 137-G34 and ESO 323-G32, both a 
TD and a RD model provide an acceptable description of the broad band spectra. An in- 
terpretation in terms of a reflection dominated continuum is preferred and considered more 
reasonable (see §4.4 and §4.5). The photon indices are in the range 1.4-1.9 and thus typi- 
cal, or possibly slightly flatter than the average values of Seyfert 1 galaxies. This suggests 
that the intrinsic continuum, when visible, is reasonably well constrained, and highlights the 
importance of the hard X-ray band in the study of heavily obscured AGN. The absorption 
corrected luminosities of the sources best fltted by the TD model are in the Seyfert range 
(-^2-10 keV ~ 10^2"^^ erg s^^). The observed hard X-ray luminosities of the reflection domi- 
nated AGN are of the order of 10^^ erg s^^. The intrinsic luminosity can be approximated 
by Lint ~ Lobs/{C X A2-10) where C = AQ/An and AQ is the solid angle illuminated by the 
central source reflecting X-rays, A2-10 is the albedo of Compton reflection in the 2-10 keV 
energy range. The albedo is a weak function of the slope of the primary spectrum and, for 
the range of slopes of the sources in our sample, is of the order of 7%, assuming reflection 
from a slab covering a solid angle of 2tt at the source. We note that the albedo due to 
torus reflection is likely to be lower than the quoted value (Murphy & Yaqoob 2009). The 
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intrinsic luminosities for the two reflection dominated Seyferts are 1.4 - 3.2 xlO^^C"^ erg 
s~^, consistent with the X-ray luminosities of Seyfert galaxies. The ratio between observed 
and intrinsic luminosity in the Circinus galaxy (Matt et al. 1999) is about a factor two lower 
than that estimated for the objects in our sample. Given that C < 1 and considering the 
uncertainties associated to the Compton reflection albedo from distant matter, the agree- 
ment is rather good. The reflection component intensity, for those sources best fitted with a 
TD model, is poorly constrained; however, it seems to be relatively weak. The iron Ka line 
equivalent widths in Compton thin sources are consistent with a transmission origin through 
the observed column densities. In the RD sources, the large EW are further supporting the 
presence of a reflection dominated continuum in these objects. An emission line feature at 
~ 7 keV (rest frame) is detected in the most obscured AGN. The best fit line energy and 
measured intensity are consistent, in all cases, with iron K/3 emission. 



4-6.2. Soft X-ray spectra 

Soft X-ray emission is ubiquitous in the sources of our sample, although the intensity 
and shape of the soft X-ray spectrum vary significantly from object to object (Fig. 1). The 
fraction of scattered flux, in the three objects whose hard X-ray continuum is best fitted 
with a TD model, is relatively low and below about 1%. These sources would be dubbed 
"hidden" AGN by Winter et al. (2008). In all sources, including the two RD objects for 
which the scattering fraction can not be defined, a soft thermal component with temperatures 
in the range 0.3-0.8 keV is also required. The 0.2-2 keV X-ray luminosity associated to the 
thermal component is in the range 1-7 xlO^'^ erg s^^, consistent with the values observed 
in "normal" galaxies (i.e. Norman et al. 2004). While the physical origin of the thermal 
component could be associated to diffuse hot plasma in the host galaxy, it is important to 
note that the best fit temperatures are derived by fitting the iron L-shell emission at 0.7-1 
keV. We also note that line-like features below ~ 2 keV are present in the residuals, with 
respect to a fit with a thermal component and a power law (Fig. 3). Prompted by these 
considerations, we explored a different possibility for the origin of soft X-ray emission. A 
photoionized plasma is found to be an excellent description of the soft X-ray spectra for a 
large sample of obscured Seyfert 2 observed with the high resolution RGS gratings on board 
XMM (Guainazzi & Bianchi 2007). The spectrum of photoionized gas is usually modeled 
with specific codes (i.e. xstar). However, the available counting statistics and the limited 
CCD energy resolution of the present data are not such to adequately constrain the physical 
properties of the emitting plasma. We adopted a simplified approach fitting the observed 
spectra with a power law plus individual narrow emission lines. The best fit line energy 
and intensities are reported in Table 4. In two sources of our sample, NGC 5728 and ESO 
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137-G34 (both Compton thick), the parameterization described above provides a formally 
better description of the soft X-ray spectra, though not statistically significant, than a power 
law plus thermal emission. A comparison of the best fit Suzaku spectra of three objects in 
the sample is shown in Fig. 3. From the analysis of the residuals, it is clear that line-like 
emission below 2 keV is not properly accounted for by a thermal plasma fit. 



5. Discussion and Conclusions 

Thanks to the broad band energy range and sensitivity of the X-ray detectors onboard 
Suzaku it has been possible to study, with good accuracy, the spectral properties of a sam- 
ple of five hard X-ray selected type 2 AGN. The presence of cold obscuring gas has been 
unambiguosly established in all the sources. Although this result was expected, given the 
sample selection criteria, the X-ray spectra show a high degree of spectral complexity. In 
particular, the absorption column density and the relative intensity of the various spectral 
components differ from object to object. 

The primary emission is seen, in three out of five sources, through column densities in 
the range ~ X0^^'^~^*^ cm~^. The fraction of the intrinsic fiux scattered into the line of sight 
and parameterized by fscatt in the partial covering fits is relatively low and the reflection 
component intensity weak. On the basis of Suzaku observations of six Swift selected AGN, 
thus an almost identical selection criterium, Eguchi et al. (2009) found that sources with 
a low scattering fraction {fscatt < 0.5 %) also have a relatively strong reflection component 
{R ^ 1). They dubbed these objects as "new" type AGN (see also Ueda et al. 2007) and 
proposed they may represent the tip of the iceberg of a hitherto uncovered population of 
AGN obscured by a geometrically and optically thick torus (Levenson et al. 2002). 

The analysis of the present data does not support the "new" type AGN interpretation. A 
possible explanation may be due to the strong degeneracy in the flt parameters. For example, 
if the reflection component is also absorbed (like in the Eguchi et al. 2009), a solution in 
terms of a stronger reflection intensity and lower normalization of the scattered component is 
found, and results similar to Eguchi et al. (2009) are recovered also in our sample (Comastri 
et al. 2009). However, we believe that the quality of our data does not allow us to constrain 
absorption on the reflection component; moreover, our best flt solution has the advantage 
to have a lower number of free parameters. It is also worth pointing out that the intensity 
of the reflection spectrum from Compton thick gas strongly depends on the geometry of 
the system. The reflection fraction, obtained from fltting disc-reflection models to Compton 
thick AGN, cannot be directly associated to a solid angle (Murphy & Yaqoob 2009). A better 
understanding of the nature of "new" type AGN and whether there are two distinct classes 
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of Seyfert 2, possibly related to the geometry of the obscuring medium, or a continuous 
distribution of spectral parameters requires a larger sample of objects with a good counting 
statistic to allow for a detailed spectral analysis. 

The observational evidence of obscured AGN with a low scattering fraction is at variance 
with the average values adopted for the above parameters in the Gilli et al. (2007) XRB 
synthesis models. We have tested the impact of different normalizations for the scattering 
component in absorbed {Nh > 10^^ cm^^) AGN. The fit to the XRB spectrum and source 
counts are basically unaffected, unless fscatt exceeds values of the order of 10 %. High fscatt 
values appears to be the exception, rather than the rule, and are definitely ruled out among 
highly obscured Swift selected AGN. 

The possibility that photoionized gas is responsible for the soft X-ray emission in obscured 
AGN is consistent with the present observations. It is interesting to note that the residuals 
in the soft X-ray band of the NGC 5728 Chandra spectrum (Zhang et al. 2006) and ESO 
137-G34 XMM-Newton spectrum (Malizia et al. 2009), both fitted with thermal emission, 
suggest the presence of unaccounted for emission lines at energies close to those reported 
here. Unfortunately, the physical status of the ionized gas cannot be constrained by the 
present observations. Only a limited number of emission lines are individually detected 
above the 3a level. It is however reassuring that the best fit energies reported in Table 4 
are consistent with those expected by ionised metals observed with the XMM-Newton and 
Chandra gratings observations of much brighter Seyfert 2 galaxies (Brinkman et al. 2002; 
Kinkhabwala et al. 2002; Guainazzi & Bianchi 2007). The ~ 0.9 keV feature, identified with 
Ne IX appears to be ubiquitous, possibly with the exception of NGC 4992 which is, however, 
the faintest soft X-ray source in the sample. 

It is worth noting that the three brightest soft X-ray sources in the sample, all of 
them Compton thick (Fig. 1), show strong [oill] ionization cones in optical images (e.g. 
Wilson et al. 1993). Moreover, extended soft X-ray emission almost spatially coincident 
with [oill] cones is detected by Chandra imaging observations of NGC 5728 (Zheng et al. 
2006). Furthermore, ESO 137-G34 and ESO 323-G32 are hosted by SO galaxies which are 
expected to have neither intense star formation activity, nor thermal X-ray emission from 
hot gas. 

Although it may be premature to discard, on the basis of the present X-ray observations, 
a thermal origin for the soft X-ray emission, we stress that the above described arguments 
would support the presence of photoionized plasma, which is also favored by high resolution 
RGS observations of bright obscured AGN (Guainazzi & Bianchi 2007). 

Long-look observations of selected sample of obscured AGN are clearly needed to elucidate 
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the nature of the soft X-ray emission in these objects, though firm conclusions can be drawn 
only by high resolution (with either gratings or calorimeters) X-ray observations. 

We thank Giovanni Miniutti for extremely useful discussions. Support from the Italian 
Space Agency (ASI) under the contracts ASI-INAF 1/088/06/0 and ASI-INAF 1/23/05 is 
acknowledged. 
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Table 1. The Sample 



Name 


z 


Exposure Time 


(ks) 


XIS" 


PIN 


NGC 4992 


0.025 


31.5 




3.30±0.06 


6.4±0.5 


NGC 5728 


0.009 


37.1 




3.10±0.06 


7.9±0.4 


ES0263-G13 


0.033 


37.4 




6.30±0.09 


4.0±0.4 


ES0137-G34 


0.009 


78.0 




1.80±0.04 


1.3±0.2 


ES0323-G32 


0.016 


80.5 




1.00±0.04 


0.7±0.2 



Note. — Background subtracted count rates in units of 10 counts 
-1 

^XIS count rates are averaged over 2 or 3 CCD units. 



O^ 



Table 2. Best fit spectral parameters 



Parameters 


NGC 4992 


NGC 5728 


ES0263-G13 


ES0137-G34 


ES0323-G32 


j^Gal (^o22 j.jn-2) 


0.02 


0.07 


0.11 


0.24 


0.08 


Nh (1022 cm-2) 


56±6 


139±15 


28±3 


1000(a) 


1000(a) 


Th 


1.56±0.13 


1.70±0.16 


1 qi+0.19 


1 43+0-24 
^•^"^-0.16 


1.62±0.45 


fscatt (%) 


<0.3 


n 7+0.5 
^■'-0.3 


<0.8 


_6 


- 


JlkeV 


<0.63 


3.3 


<2.9 


4.3 


1.9 


R 


3+"-2 


0.2±0.1 


5+"-^ 
U.0_o4 


OO (a) 


oo(a) 


E{Ka) (keV) 


f, OO+0.02 


6.39±0.01 


6.39±0.04 


6.40±0.01 


6.42±0.02 


EW (kcV) 


0.34±0.07 


1 no+0.39 
^■'J'^-0.40 


0.08±0.04 


i-46t?,:^^ 


1.93±0.28 


E{Kl3) (keV) 


- 


7.05±0.04 


- 


6.98J:H? 


6 97+"-"» 


EW* (keV) 


- 


n7+""5 


- 


f,o,+0.12 


0.26±0.15 


kT (keV) 


34+0-35 


0.28±0.03 


n co-i-0.08 

U.0^_0 24 


0.84±0.06 


0.78±0.20 


F (0.5-2 keV) 


0.02 


0.16 


0.07 


0.14 


0.05 


F (2-10 keV) 


2.4 


1.5 


3.8 


0.53 


0.40 


L (2-10 kcV) 


1.6 


0.48 


2.7 


0.01(c) 


0.023(e) 


xVd.o.f. 


126.0/123 


169.5/143 


177.8/207 


169.7/139 


79.0/89 



~J 



erg s . 

''For the sources best fitted with a RD model, the column density is fixed to 10^^ cm~^and the ratio 
of the reflected to the transmitted component diverges by definition. 

''No acceptable solutions arc found linking the two power law indices. Leaving both the power law 
slopes free to vary, the best fit in the soft X-rays is 2.45±0.40. 

'^Fluxes at 1 kcV in units of lO^'"* ph cm^^ s^^ kcV^"'. 

"^Thc lines EW is computed against the total continuum at the best fit line energy. For the sources best 
fitted by a TD model, the line is also absorbed by the quoted column density. 

''Observed luminosities. 



Table 3. Best fit spectral parameters 



Parameters 


NGC 4992 


NGC 5728 


ES0263-G13 


ES0137-G34 


ES0323-G32 


j^Gal (1022 cm-2) 


0.02 


0.07 


0.11 


0.24 


0.08 


Nh (1022 cm-2) 


55±6 


140t| 


26±1 


1000(a) 


1000(a) 


Th 


i.t>D_(, ;^4 


1 fi«+0.05 
-'-■O»_o.04 


1 70+0.05 

^-■'^-0.04 


1 r;s+0.20 
i.0»_0.16 


1 oc+0.47 
i.o0_0.40 


n 


...(c) 


2.7±0.3 


2.40±0.90 


n 11+0.28 

•^■^^-0.24 


...(d) 


JlkeV 


<0.83 


6.3 


3.1 


12.6 


2.3 


R 


0.32±0.12 


n •54+0.04 


< 0.68 


00 (a) 


00(a) 


E(Ka) (keV) 


f. 00+0.02 


6.39±0.01 


6.39±0.05 


6.40±0.01 


6.42±0.02 


E.W. (keV) 


0.35±0.07 


1 90+0.14 
l-28_o.i3 


n ns+0.045 
0-08-0.035 


1.6±0.2 


2.0±0.3 


E{KI3) (keV) 


- 


7.05±0.05 


- 


fi qq+0.06 


6.97±0.08 


E.W. (keV) 


- 




- 


•H^+0.12 


0.27±0.17 


x'/d.o.f. 


120.3/119 


135.3/133 


172.7/206 


146.7/131 


78.8/87 



Note. — Spectral fits parameters obtained when the soft X-ray spectrum is fitted with a power law 
plus individual Gaussian lines. 

''For the sources best fitted with a RD model, the column density is fixed to 102^ cm^2and the ratio 
of the reflected to the transmitted component diverges by definition. 

''Photon index of the power law representing the photoionized continuum. 

'^The continuum flux is extremely weak and formally consistent with zero. 

''The power law slope in the soft band is unconstrained and fixed to be equal to the hard power law 
slope. 

''Fluxes at 1 kcV in units of 10"'°' ph cm~2 s~' keV~'. 



00 



Table 4. Soft X-ray lines 



Parameters 


NGC 4992 


NGC 5728 


ES0263-G13 


ES0137-G34 


ES0323-G32 




Line ID" 


El (keV) 


U-'U_o.02 


0.65±0.02 


0.69±0.03 






VIII 


(0.65) - Fe xvii (0.72) 


Fluxi 


0.9±0.5 


3.7±1.0 


1 9+2-2 










E2 (keV) 


0.84±0.02 


O.SSiO.Ol 


0.87±0.02 


0.91±0.02 


0.94±0.02 


Ne IX (0 


90-0.92) - Fe xviii (0.86) 


Flux2 


0.4±0.2 


1.0±0.2 


0.94±0.3 


0.67±0.30 


0.3±0.1 






E3 (keV) 




0.98±0.02 




1.03±0.01 






Ne X (1.02) 


FIUX3 




0.4±0.1 




0.90±0.20 








Ei (keV) 


1 fi9+0-06 
^■''^-0.03 


1.76±0.02 




1.79±0.02 


1.81±0.03 




Si XIII (1.84) 


Flux4 


0.09±0.04 


0.11±0.04 




0.18±0.05 


0.08±0.04 






Ei (keV) 




2.35±0.03 




2.53±0.05 






S XV (2.44) 


Fluxs 




0.15±0.06 




n 1^+0.03 








Eq (keV) 




3.13±0.04 




3.11±0.02 






Ar xvii (3.11) 


Fluxe 




O.llitO.05 




0.10±0.04 









Note. — Emission-line fluxes in units of 10 



Errors on both line energy and flux are at Icr. 



CO 



^Suggested line identification and corresponding energies in keV (from Guainazzi & Bianchi 2007 and House 1969). The approximate 
energy interval for the Ne ix triplet is reported. 
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Fig. 1. — The broad band Suzaku spectra of the five sources, covering about two decades in 
frequency, in a vYy representation. The Y-axis normahzation is arbitrary. 
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Fig. 2. — The Suzaku broad band spectra and data to model ratio for the 5 Seyfert galaxies 
of the sample. In the left column we plot the sources best fitted with a TD model, and in 
the right column the two reflection dominated X-ray spectra. The spectral parameters of 
the broad band fits shown in the figures are reported in Table 2. 
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Fig. 3. — The Suzaku soft X-ray spectra (< 4 keV) and data to model ratio of NGC 5728, 
ES0137-G34 and NGC 4992 (from top to bottom). In the left column the best fit to the 
soft X-ray spectrum obtained with a power law plus a thermal model is reported. In the 
right column the best fit with a power law plus Gaussian lines is shown. For NGC 4992, 
only the soft X-ray spectrum, without residuals, is reported. 



